Particles and Anti-Particles in a Relativistic Bose Condensate 



Luca Salasnich 

Istituto Nazionale per la Fisica della Materia, Unitd di Milano Universitd, 
Dipartimento di Fisica, Universitd di Milano, 
Via Celoria 16, 20133 Milano, Italy 
E-mail: salasnich@mi.infm.it 

Abstract 

We study the Bose-Einstein condensation (BEC) for a relativistic ideal gas 

of bosons. In the framework of canonical thermal field theory, we analyze 

the role of particles and anti-particles in the determination of BEC transition 

temperature. At the BEC transition point we obtain two universal curves, i.e. 

valid for any mass value: the scaled critical temperature as a function of the 

scaled charge density of the Bose system, and the density ratio of anti-particles 

versus the scaled critical temperature. Moreover, we numerically investigate 

charge densities and condensed fraction ranging from the non-relativistic to 

the ultra-relativistic temperature, where analytical results are obtained. 
PACS: 05.30.Jp, ll.30.Qc 
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I. INTRODUCTION 



Nowadays more than twenty experimental groups have achieved Bose-Einstein conden- 
sation (BEC) in clouds of confined alkali- metal atoms. ^ These studies have renewed the the- 
oretical interest on non-relativistic^ but also relativistic BEC. In the past years, relativistic 
BEC has been analytically investigated by several authors using Euclidean-time functional 
integration^"^ but a quantitative numerical analysis with the temperature ranging form the 
non-relativistic to the ultra-relativistic limit has never been performed. 

In this paper we consider BEC in the case of a relativistic non-interacting Bose gas de- 
scribed by a complex scalar field (f){x). We derive the exact equation of motion of in 
the grand canonical ensemble of equilibrium statistical mechanics® without invoking func- 
tional integration^ but using instead canonical field theory. By means of the Bogoliubov 
prescription^'^, we write down the equations of the condensate order parameter and of ther- 
mal particles and anti-particles. We numerically investigate the effect of particles and anti- 
particles in the determination of BEC transition temperature Tc- Moreover we study the frac- 
tion of anti-particles in the system as a function of temperature and derive ultra-relativistic 
formulas for Tc and the condensed fraction in a generic d-dimensional space. 

II. SCALAR FIELD AND LEGENDRE ANTI-TRANSFORMATION 

The Lagrangian density of a non-interacting complex scalar field (f){x) is given by 



where m is the mass of the identical bosons described by the scalar field. To study the 
finite-temperature properties of a field-theory one needs the Hamiltonian H of the system. 
In our case, the canonical conjugate momentum Il{x) of the scalar field is 




(1) 




,+ 



(2) 



and the Hamiltonian density reads 
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n = n+n + V0+ • V0 + m^0+(/) . (3) 

Note that the invariance of the system under a global gauge U(l) field transformation implies 
the conservation of the electric current and the conserved charge density is 

Q = i (0+0 - 00+) = i (n+0+ - n0) . (4) 

Equilibrium statistical mechanics tells us that the grand canonical partition function Z 
of a quantum system of Hamiltonian H and conserved charge Q is given by 



Z^Tr 



^-{H-nQ)/T 



(5) 



where T is the temperature of the thermal reservoir and /j, is the chemical potential.^ Thus 
the system in the grand canonical ensemble is described by an effective Hamiltonian density 
Tin — H — /iQ, namely 

= n+n + V0+ • V0 + y (0) - i//(n+0+ - n0) . (6) 

To find the effective Lagrangian density we observe that the Hamilton equation for is given 

by 

= ^ = n+ + i/.0. (7) 

Using the Legendre anti-transformation one can easily obtain the effective Lagrangian den- 
sity Cn = 0+n+ -|- 0n — Tii^, namely 



Cl,^jC + tfM (0+0 - 00+J + /0+0 . (8) 

Therefore, the introduction of the chemical potential /i is equivalent to the use of an effective 
Lagrangian that can be obtained with a shift 

d d . 

in the time partial derivative of the bare Lagrangian C. It is important to observe that the 
same result could be obtained by means of the Euclidean-time functional integration^ and 



that the shift found holds also in the fermionic case.^^ Finally, the Euler-Lagrange equation 
of the effective Lagrangian density reads 



0, 



(10) 



at 



where □ = ^ — is the d'Alambert operator. 



III. BOSE-EINSTEIN CONDENSATION AND BOGOLIUBOV PRESCRIPTION 



In a Bosonic system one can separate Bose-condensed particles from non-condensed ones 
by means of the Bogoliubov prescription^^ that is given by 



= $ + 77 , 



(11) 



where 



(12) 



is the order parameter of the Bose condensate (a classical complex scalar field), namely 
the non- vanishing thermal average of the Bosonic field, and r]{x) is the operator of the 
non-condensed or thermal particles, such that (77) = 0. 

The exact equation of motion of the order parameter $(x) is obtained by calculating the 
thermal average over the equation of motion of the scalar field 4>{x). If we have a static and 
homogeneous order parameter, then 



(13) 



from which it follows that there is macroscopic occupation of the lowest single-particle state 
($ ^ 0) only if |//| = m. 

The exact equation of motion of the fluctuation operator r]{x) is easily obtained by 
subtracting the exact equation of ^{x) to the equation of (j){x). Note that, in the non- 
interacting case, the equations of (f){x), ^{x) and r]{x) are formally identical but they have 
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different meanings. In addition, by using the Bogoliubov prescription and the effective 
Lagrangian, the grand canonical charge density of the system reads 

Q = Qo + Q , (14) 

where 

Qo = 2//|$|2 (15) 
is the contribution due to the Bose condensate and 

Q — i (rj'^r] — rjr)'^^ (16) 
is the contribution due to thermal particles and anti-particles (see also Ref. 3). 

IV. BOSE CONDENSATE AND THERMAL PARTICLES 

The relativistic complex scalar field operator ri{x) satisfies the equal-time commutation 
rule 

[r]{^,t),r)+{y,t)]^S'{^-y). (17) 

The operator r){x) can be Fourier decomposed into a single-particle basis of particles and 
anti-particles 

( j(k-x-t<jfct) -i(k-x-(Iifet) \ 

-7^^^^+ bt] , (18) 

where we have used the symbols cuk and Cuk to indicate the value of energy for particles and 
anti-particles, respectively. The Bose operators for particles and anti-particles satisfy the 
canonical commutation relations 

[ok, a^,] = [6k, b^,] = 5kk' , (19) 

and all other commutators are zero. In addition, one imposes the following Bose-Einstein 
thermal averages 



The energies uik and cDfc are determined by inserting the decomposition of the field rj{x) in 
its equation of motion, namely Eq. (10) with r]{x) instead of (t){x). In this way, one finds 
two decoupled algebric equations: 

ujI - 2nujk + pL^ - - k"^ ^ , (21) 

Ol + 2iiujk + -m'^ - k'^ = Q , (22) 
which give the physical solutions 



ijjk = vk^ + m? — n , oJk = vW+m? + ji . (23) 

The Fourier decomposition and the energies ui^ and uj^ enable us to calculate the thermal 
average q— (Q) of the non-condensed charged density, which is given by 

^ ~ ^ p(Vfe^+m^-M)/T _ 1 ~ p(Vfc^+m^+M)/T _ 1 ' ^^4) 

Thus, q — rii — n2, where rii — Z^k(ctk^k) is the average density of particles and n2 — 
Z]k(^k^k) is the average density of anti-particles. Note that 77-2 is obtained from rii with 
the substitution — > — /x. The chemical potential /j, describes both bosons and antibosons: 
the sign of /i indicates whether particles outnumber antiparticles or vice versa. Moreover, 
because both rii and 77-2 must be positive definite, it follows that < m. 

Obviously, the total number of particles is not conserved because of the production of 
antiparticles, which becomes relevant when T is comparable with m. The conserved quantity 
is the net charge density q — qo + Q, where q'o = Qo = 2/x|0p is the condensed charge density 
and q = (Q) = rii —n2 is the the difference between the density of particles and the density 
of anti-particles. The condensed charge density qo is non-zero only below the BEG transition 
temperature Tp. The condensed charge density corresponds to A; = in Eq. (24) and it is 
thus given by 
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where and n2 are the density of condensed particles and anti-particles, respectively. It 
is easy to show that 

hmnf\T) = qo, lim4°)(T) = . (26) 

As expected, at T = the particles are all in the condensate and there are no anti-particles. 
Moreover, in the hmit T — > the asymptotic behavior of the chemical potential and of 
the density 712'^ of condensed anti-particles read 

Mr)~m-ri,i('*-ti'j , 

To our knowledge, this is the first paper where these simple asymptotic relations have been 
explicitely written down. 



V. NUMERICAL AND ANALYTICAL RESULTS 

The behavior of the Bose gas ranging from the non-relativistic to the ultra-relativistic 
regime can be numerically investigated by means of the non-condensed thermal charge den- 
sity q given by Eq. (24). In particular, we work in the thermodynamic limit substituting 
the sums in Eq. (24) with integrals and find 



1 

dk 



27r2 







1 1 



g(\/fc^+m2-/i)/T _ I g(Vfe^+m2+/i)/T _ 



1 



27r2 







dk [ni{k) — n2{k)] , 
(28) 



where ni{k) and 712 (/c) are the density profiles in momentum space for particles and antipar- 
ticles, respectively. It is important to stress that the previous formula has an useful scaling 
property: the chemical potential, the temperature and the momentum can be measured in 
units of m and the densities in units of m^. This follows from the fact that, with h — c— 1, 



the mass, the chemical potential, the temperature and the momentum have the same unit: 
the energy, while the length is measured in units of the inverse of energy. 

We first consider the case T > Tc. Given the charge density and the temperature, the 
chemical potential is fixed by Eq. (28). In such a way one determines also the fraction of 
anti-particles in the system. As shown in Fig. 1, where we plot 77-2/711 as a function of the 
scaled temperature T/m for different values of the scaled charge density q/m^ of the Bose 
system, one can identify two regimes: the non-relativistic regime {T/m <^ 1) and the ultra- 
relativistic regime {T/m ^ 1). In the non-relativistic regime, the fraction of anti-particles of 
the system is negligible. In the ultra-relativistic regime the fraction of anti-particles becomes 
relevant. The curves with fixed charge density q — q {qo — 0) ends at the scaled critical 
temperature T^/m, where the Bose- Einstein condensate appears. The role of temperature 
in the formation of anti-particles is also shown in Fig. 2, where we plot the density profiles 
in momentum space of particles and anti-particles for three increasing values of the scaled 
temperature. 

We observe that in the ultra-relativistic regime one can derive analytical results by 
performing a Taylor expansion of q at first order in /i. After straightforward but tedious 
calculations one finds 

r^l = #T^ + ^T^ m = ^-^T' - f^T\ q^^T\ (29) 

TT^ 6 TT"' 6 3 

where C,{x) is the Riemann ^-function. These analytical results, confirmed by our numerical 
calculations, show that although n2/n\ — > 1 as T — >^ cxo, the charge density q = n\ — n2 goes 
to infinity. 

We have previously shown that the critical temperature Tc at which BEG occurs corre- 
sponds to = m. At the BEG transition temperature the thermal charged density q 
can still be determined from Eq. (28). In fact, by inverting the function q{Tc,m — /i) one 
finds the transition temperature. In Fig. 3 we plot two curves which do not depend on the 
value of the mass of Bosons in the gas (we call them universal curve). The first unversal 
curve is the scaled critical temperature Tc/m as a function of the scaled charge density 
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q/rri^. The second universal curve is the ratio 712/711 between anti-particles and particles 
as a function of the scaled critical temperature Tc/m. As expected, grows with q and 
n2/ni. Moreover, at a fixed density ratio n2/ni, it is easier to get high-temperature BEC 
with heavy-mass particles. The two universal curves of Fig. 3 can be compared with ultra- 
relativistic analytical results. Prom (29) one immediately finds that the critical temperature 
(/X = m) is given by 

and the density ratio reads 



^0= - , (30) 
\mj 



„ C(3)t3 _ £^T2 

- = ci.^ ' ■ (31) 

^ 7r2 -'cTg-'c 

Fig. 3 shows that while the formula of the density ratio n2/ni is valid only in the ultra- 
relativistic region (for Tc ^ it predicts the wrong limit 712/712 —1), the formula of the 
critical temperature is quite accurate also at low temperatures. Note that the formula of 
the critical temperature has been first obtained by Kapusta.^ We now extend it to the case 
of a ultra-relativistic gas in d-dimensional space. 

The charge density our system of non-interacting bosons can be re-written as 

/■oo V \ \ -\ 

q = de p{e) [ ^(,_,)/^ _ ^ - ^(e+,)/T _ J ' (32) 
where p(e) is the density of states. It can be obtained from the formula 

P(e) = J-0j-A'-H{k)), (33) 
where H{k) is the classical single-particle Hamiltonian of the system in a d-dimensional 



space. The classical single-particle Hamiltonian of a relativistic ideal gas is H — ^/k^ + m? 
and the density of states reads 

In the ultra-relativistic limit 

the density of states is simply p(e) = {2n'^/'^)/{{2nYT{d/2))e'-'^-'^\ In this case, by using 
again the Taylor expansion of q at first order in /i with T — Tc one finally obtains 
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( 



(2nYr{d/2) q 
47r<^/2r((i)C(ci - 1) m 



) 



i/(d-i) 



(35) 



where r{x) is the factorial function. Because C(l) = oo, it follows that for a homogeneous 
relativistic gas there is BEC only for c? > 2, as in the case of a non-relativistic homogeneous 
gas (see also Ref. 12, Ref. 13). 

Below Tc, a macroscopic number of particle occupies the single-particle ground-state of 
the system (go 7^ 0). The Eq. (28) gives the charge density q — q — qq of non-condensed 
particles. In this way, form Eq. (15), one determines the order parameter that is such 
that 



Thus, the condensed fraction qo/q can be numerically calculated as 1 — q{T)/q{Tc). In Fig. 
4 we show the condensed fraction as a function of the temperature for different values of the 
scaled charge density q/m^ of the Bosonic particles. In the ultra-relativistic regime, from 
(28) and (29) one finds 



namely the condensed fraction has an inverted-parabola shape. It is important to stress 
that, although we are able to determine the charge density qo of the Bose condensate, our 
formalism cannot tell us the fraction of anti-particles into the condensate. 

For an ideal gas of charged massless Bosons (m = 0) it follows, form Eq. (30) and 
Eq. (35), that Tc — 00 and go = Q' at any temperature, all net charge resides in the Bose 
condensate. Nevertheless, if the thermal average of the charge is not conserved (// = 0), 
i.e. a gas of photons, then BEC does not take place. Finally, by using the previously 
discussed procedure, one finds that the Bose condensed fraction for a ultra-relativistic gas 
in a d-dimensional space is given by 



g(T,) - q{T) 



(36) 



2m 




(37) 




(38) 



remembering that Tc ^ 00 as 0? ^ 2. 



10 



VI. CONCLUSIONS 



We have studied thermal properties of a non-interacting relativistic Bose gas by analyzing 
in detail the fraction of anti-particles in the system. By using a finite-temperature operator 
formalism, we have obtained the equation of the Bose condensate, described by a complex 
classical order parameter, and the equation of non-condensed particles and anti-particles. At 
zero temperature the particles are all in the condensate and there are no anti-particles. In 
the limit of zero temperature we have determined the asymptotic behavior of the density of 
condensed anti-particles. The charge density and the density of particles and anti-particles 
have been analyzed as a function of temperature raging from the non-relativistic to the ultra- 
relativistic regime. We have determined two universal curves at the BEC transition point: 
the scaled critical temperature Tc/m as a function of the scaled charge density q/rm? of the 
Bose gas, and the ratio between anti-particles and particles as a function of the scaled critical 
temperature Tc/m. Moreover, we have investigated the condensed fraction as a function of 
the scaled temperature for increasing values of the scaled charge density of the gas. Finally, 
analytical results have been found in the ultra-relativistic region. Our analytical formulas for 
a ultra-relativistic Bose gas in d-dimensional space generalize previous findings with d = 2>. 

In conclusion, we observe that detailed analytical and numerical investigations can be 
also performed in the case of an interacting relativistic Bose gas, at least in the Bogoliubov- 
Popov mean-field approximation. This is one of our future projects. 
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FIG. 1. Density ratio n2/ni between particles and anti-particles vs scaled temperature T/m. 
Curves for different values of the scaled charge density q/m^ of the Bose gas above the critical 
teperature Tc, where each curve ends. 
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FIG. 2. Density profile in momentum space for particles (full line) and anti-particles (dashed 
line). Scaled charge density of the Bose gas: qjvp? = 0.1. Prom left to right: T/m = 1, T/m = 1.5, 
T/m = 2. 
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FIG. 3. Universal curves at the BEC transition point. Scaled critical temperature Tc/m vs 
scaled charge density q/m^ (above). Density ratio 712/ n\ vs scaled critical temperature Tc/m (be- 
low) . Full lines are numerical results and dashed lines are analytical results in the ultra-relativistic 
limit. 
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FIG. 4. Condensed fraction q^/q vs scaled temperature T/m. Curves for different values of 
the scaled charge density q/m^ of the Bose gas below the critical temperature. 
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